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The magnetotransport properties of antidot lattices containing artificially designed grain bound- 
aries have been measured. We find that the grain boundaries broaden the commensurability res- 
onances and displace them anisotropically. These phenomena are unexpectedly weak but differ 
characteristically from isotropic, Gaussian disorder in the antidot positions. The observations are 
interpreted in terms of semiclassical trajectories which tend to localize along the grain boundaries 
within certain magnetic field intervals. Furthermore, our results indicate how the transport through 
superlattices generated by self-organizing templates may get influenced by grain boundaries. 

PACS numbers: 73.23.-b, 73.63.-b 
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I. INTRODUCTION 

Superlattices imposed on two-dimensional electron 
gases (2DEGs) by lateral patterning have revealed many 
fascinating effects over the past two decades. One version 
of such superlattices are antidot lattices, where the po- 
tential maxima lie above the Fermi energy fii^i^ The longi- 
tudinal magnetoresistance shows resonances at magnetic 
fields where the cyclotron orbit is commensurate with 
the lattice^, an effect which originates from the classi- 
cally chaotic character of such systems, as demonstrated 
by simulations based on the Kubo formalism,^'^ This 
model has also been used to explain the nonlinear Hall 
resistance of antidot lattices)^ which can even develop a 
pronounced negative Hall resistance at small magnetic 
fields.^ Many quantum effects in lateral superlattices 
have been reported as well. For example, i?-periodic 
oscillations, superimposed on the commensurability os- 
cillations, have been explained by quantization along 
closed orbits^'* in square lattices. In hexagonal lattices, 
a transition from Altshuler-Aronov-Spivak oscillations at 
small magnetic fields to Aharonov-Bohm oscillations at 
larger fields has been observedi^iiSiii Moreover, the long- 
standing prediction of the Hofstadter butterfly)^ the 
fractal energy spectrum of a periodic potential in strong 
magnetic fields, has been experimentally confirmed in 
weak periodic superlattices J^ii^ii^ Antidots have further- 
more been used to demonstrate the quasi-particle char- 
acter of composite fermions in a very intuitive wayJ^ 

These results have generated the quest to drive lat- 
eral superlattices deeper into the quantum regime which 
requires smaller lattice constants. With established 
techniques reaching their limits around periods of 100 
nmfi^ii^ self-organizational schemes offer themselves as 
an alternative. Potential techniques include masks from 
porous aluminEfcii or from diblock copolymersi^i. Re- 
cently, Melinte et al. demonstrated that a diblock copoly- 
mer system in the spherical phase can be used to gener- 
ate a two-dimensional lateral superlattice in a Ga[Al]As 
heterostructure with lattice constants below 40 nm^J^ 



Diblock copolymer systems in the hexagonal columnar 
phase have already been used to pattern, e.g., ferro- 
magnet/normal metal hybrid systemgSS as well as semi- 
conductor surfaces. ^^ However, in many self-organizing 
systems, domain formation and the corresponding disor- 
der at the grain boundaries (GBs) poses a limitation to 
the periodicity of the patterns. ^^'^^ Based on this back- 
ground, the question how GBs in two-dimensional super- 
lattices influence the transport properties is a relevant 
one which has not yet been addressed, to the best of our 
knowledge. Beyond that, the study of such mesoscopic 
GBs is of fundamental relevance. First of all, GBs in 
natural polycrystalline conductors usually trap impuri- 
ties, which may be charged and dominate the influence 
of the GBs on the conduction electrons i22i2^ The artifl- 
cially designed, mesoscopic GBs in lateral superlattices, 
on the other hand, are free of such impurities and there- 
fore allow us to study the purely geometrical impact of 
a GB on the transport properties. Second, the electronic 
wave packets in antidot lattices can be highly localized 
on the scale of the lattice constant and thus within the 
GB, a situation which is not possible in conventional con- 
ductors. 

Here, we report an investigation of artificial but real- 
istic GBs defined in a hexagonal antidot lattice. The 
GBs have been generated by Monte-Carlo simulations 
with appropriate boundary conditions, followed by a pat- 
tern transfer into the 2DEG via conventional lithographic 
techniques. We observe a surprisingly weak, but never- 
theless characteristic, anisotropic damping and shift of 
commensurability resonances in the magnetoresistance. 
These experimental results are supported by numerical 
simulations within the Kubo formalism. The calculations 
reveal that the resistance anisotropy has their origin in a 
small, magnetic field dependent fraction of electron tra- 
jectories that tend to localize along the GBs over large 
magnetic field intervals. 

The paper is organized as follows. In Section II, we 
present the sample preparation and the experimental 
setup. Section HI contains the experimental results and 



their interpretation, 
conclusions. 



Section IV gives a summary and 
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II. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

Conventional modulation-doped GaAs/Alo.2Gao.8As 
heterojunctions with a two-dimensional electron gas 
(2DEG) 100 nm below the surface were grown by molec- 
ular beam epitaxy. The pristine 2DEG has an elec- 
tron density oi n — 3.2 x 10^^ m'~^ and a mobility of 
/x = 68 m^/Vs, corresponding to an elastic mean free path 
of I = 6.4 /im and a Drude scattering time of td = 26 ps 
at liquid helium temperatures. Since all components of 
the (anisotropic) resistivity tensor must be accessible, we 
defined a square (20 /im x 20 /im) van der Pauw-geometry 
with 8 contacts by electron beam lithography and sub- 
sequent ion beam etching with low-energy (250 eV) Ar"*" 
ions.S4 The etch depth was 20 nm which is sufficient to 
deplete the 2DEG. For Ohmic contacts to the 2DEG we 
alloyed Ni/AuGe pads in the heterostructure. 

The 2D lattice containing the GB was numerically gen- 
erated by using a standard Metropolis Monte-Carlo sim- 
ulation algorithm with periodic boundary conditions and 
a repulsive, two-dimensional Yukawa-Potential which is 
cut off after three lattice constants. ^^ Without any initial 
conditions (i.e. seeds), this would generate a perfectly 
hexagonal lattice. In order to enforce the formation of a 
GB, we defined two seeds, each consisting of 19 sites ar- 
ranged in perfectly hexagonal geometry, separated by 20 
lattice constants and rotated to each other by an angle 
oi 9 — 21.8 °, see Fig. HJa), which corresponds to the E7 
coincidence site lattice within one of the standard classi- 
fication schemes for GBa^^i^l, meaning that one out of 7 
lattice sites coincide in an overlay of the two grains. It be- 
longs to the set of energetically favorably GBs according 
to various modelSf^S and occurs frequently in real sys- 
tems see, e.g. Puglisi et al.'^^. Fig. 3(a) or Black et alf?°. 
Fig. 1(c). We believe, however, that the experimental 
results reported below are rather insensitive to the kind 
of coincidence site lattice, and we note that we are un- 
aware of a scheme for classifying GBs according to the 
microscopic details of the local geometry at the bound- 
ary, which would probably be most appropriate for our 
experiments. 

The simulated lattices with GBs were transferred to 
the pattern generator of an electron beam writer. The 
GBs were oriented parallel to the mesa edge in y- 
direction. The superlattices were patterned by electron 
beam lithography and ion beam etching as described 
above. The average lattice constant was a = 600 nm 
with a lithographic antidot diameter of duth ~ 200 nm. 
We opted for patterning antidot lattices with a large ra- 
tio deff/a « 0.7, where deff is the effective electronic 
antidot diameter, in order to amplify the effects due to 
the disorder at the GBs as compared to lattices with 
smaller d^ff/a ratios, at the expense of the number of 
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FIG. 1: (a) A section of the hexagonal lattice generated by 
the Monte-Carlo simulation containing a E7 grain boundary, 
characterized by the angle 6 = 21.8°. (b) Top view of the 
patterned sample prior to the deposition of the top gate, as 
seen in a scanning electron microscope. A stripe of the lattice 
around the GB along the y direction is replicated across the 
van der Pauw square in the center of the mesa that contains 
the 2DEG (light gray), (c) Histogram of the nearest-neighbor 
distances of the designed GB lattice, measured in units of the 
average lattice constant a. For this purpose, the replication of 
the GB lattice has been extended to 10* antidots. A Gaussian 
fit gives a mean value of 1 with amplitude and o" = 0.10 in 
units of a as fit parameters. 



observable commensurability oscillations in similarly pat- 
terned arrays.-il Finally, on two of the samples, a Ti/Au 
(20 nm/200 nm) top gate was evaporated. 

The magnetoresistance was measured in a helium gas 
flow cryostat with a variable temperature insert and a 
superconducting magnet with a maximum field of _B = 
8 T which is applied perpendicular to the plane of the 
2DEG. For temperatures at 100 mK, a ^He/^He dilution 
refrigerator was used. All measurements were performed 
in a four-probe setup using standard lock-in techniques 
with a current of 100 nA at frequencies of 13.7 Hz. 



We did not observe any GB effect on antidot lattices 
containing just one GB along the y-direction as indicated 
in Fig. [UJa) (not shown). Therefore, the simulated GB 
was periodically replicated without disturbing the reg- 
ular hexagonal lattice parts outside the grain boundary 
area, see Fig. [Hb). The average distance between two 
GBs was dcB ~ 2.5 /im. This is about 40% of the mean 
free path of the unpatterned 2DEG and well above the 
effective mean free path in the antidot lattice, estimated 
from the resistivity at zero magnetic field to w 280 nm . 
Therefore, we can exclude the possibility of commensura- 
bility effects due to the one-dimensional GB superlattice. 
At the same time, dcB is larger than the spatial exten- 
sion of a GB, such that its character is retained. Fig. 
[TJc) shows the distribution of nearest neighbor (NN) dis- 
tances in the GB lattice. The NN distances show marked 
deviations from Gaussian disorder. Compared to a Gaus- 
sian NN distribution with the same average distance a 
and the same single standard deviation of 0.1a, a long 
tail of large distances is observed, while the majority of 
the NN distances is slightly reduced, resulting in a his- 
togram peak at a NN distance of w 0.96a, in combination 
with a suppression of very small NN distances. This in- 
dicates that within the domain wall, a few antidots with 
large NN distances are present wich compress the sur- 
rounding lattice locally. Very small NN distances below 
K, 0.9a, however, are suppressed, due to the repulsive po- 
tential used in the Monte Carlo simulation. As a result, 
the maximum of the measured NN distribution is larger 
than the corresponding Gaussian distribution, and the 
peak appears significantly sharper around its maximum. 

Four samples of this geometry and of nominally iden- 
tical values diith/a have been measured. They all show 
qualitatively very similar behavior. Here, we show data 
from the two samples where the effects were most pro- 
nounced. In addition, two single domain hexagonal lat- 
tices in the same sample have been measured for control 
purposes. 



III. EXPERIMENTAL RESULTS AND 
INTERPRETATION 

In Fig. m^a), the longitudinal magnetoresistances Rxx 
and Ryy for a single domain hexagonal antidot lattice and 
for the lattice containing the GBs (as reproduced in Fig. 
[1]) are shown. For better comparison all measurements 
were normalized to their zero-field value. 

The measurements show the characteristic commensu- 
rability oscillations of antidot arrays. Within the sim- 
ple commensurability picture, the cyclotron radius of the 
electrons roughly matches the orbital radius around 1 and 
7 antidots at B « 230 mT and « 100 mT, respectively.^^ 
The resonance around 3 antidots is weakly pronounced 
and hidden in the wing towards lower magnetic fields of 
the resonance around one antidot. The suppression of 
this resonance is also observed in our simulations for an- 
tidots with large d/a ratios and is a simple geometrical 



effect. We also observe a weak resonance at i? w 450 mT 
which corresponds to electrons caught in between three 
adjacent antidots, as possible in arrays with large d/a 
ratios:^ Moreover, the characteristic overall resistance 
decrease as a function of increasing B is observed, which 
ends roughly where the cyclotron radius equals the spa- 
tial extension of the electron gas between two adjacent 
antidots. This was used to estimate the effective antidot 
diameter to dgff ~ 460 nm up to de// ~ 520 nm depend- 
ing on the sample and the top gate voltage. 
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FIG. 2; (a) Upper part: longitudinal magnetoresistances of a 
single domain hexagonal antidot array and a GB superlattice. 
Both samples were ungated. Lower part: the corresponding 
resistance anisotropy Ap as defined by Eq. [T] (b) Tempera- 
ture dependence of Ap for a second sample containing a GB 
lattice, (c) Comparison of the experimentally determined Ap 
of the GB lattice shown in (a) with the simulations. 

The observation of the commensurability oscillations 
demonstrates their resilience with respect to GBs. How- 
ever, compared with the single domain lattice, the com- 
mensurability peaks in the GB lattice are significantly 
damped and somewhat broadened. The peak position 



in the GB lattice that corresponds to the cyclotron mo- 
tion around a single antidot is shifted to larger magnetic 
fields, while the positions of the other two resonances are 
approximately constant. This peak shift is about 5% for 
Ryy and 10% for Rxx, and cannot be explained with de- 
viations of the lattice constants from the designed values, 
since the accuracy of the electron beam writer allows for 
1% deviation at most. Rather, we interpret this displace- 
ment as a consequence of the NN distribution. The sharp 
peak of NN distances at a reduced value corresponds to a 
locally reduced lattice constant, which is reflected in the 
position of the resonance around a single antidot while 
the resonance around 7 antidots is less sensitive to this 
effect since it averages over more NN distances. 

The single domain lattice shows isotropic behavior to 
a good approximation, whereas in the lattice contain- 
ing the GBs, the transport is anisotropic. All three 
resonances are slightly suppressed in the y-direction, 
i.e., along the grain boundaries, as compared to the x- 
direction (perpendicular to the GBs). Especially the the 
commensurability resonance at i? w 240 mT is shifted 
to smaller magnetic fields in the y direction. Further- 
more, Ryy is significantly lower than R^x at and above 
the resonance around a single antidot. We emphasize 
that this phenomenology is observed in all our GB lat- 
tices. This anisotropy somewhat resembles that one re- 
ported by Tsukagoshi et al.^'^, who measured rectangular 
antidot lattices with uni-directional, Gaussian disorder 
and observed a damping of the commensurability oscilla- 
tions as well as a shift to smaller magnetic fields for the 
current perpendicular to the disorder direction. In anal- 
ogy to this behavior, one is tempted to conclude that in 
the GB lattice, the disorder perpendicular to the GBs 
is greater than parallel to them. However, we could not 
think of a way to check this, e.g. by looking at the x- and 
y-components of the NN distributions, since our single 
domain stripes have different orientations with respect 
to each other and many antidots at the grain boundaries 
cannot be unambiguously attributed to one domain. 

We define the resistivity anisotropy Ap as 
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In the van der Pauw geometry used, Ap is equivalent to 
a correspondingly defined resistance anisotrop y'^^i'^^ and 
nicely illustrates the anisotropy effects described above. 
In the lower part of Fig. [^a) and in Fig. ^h), Ap{B) 
for our measurements is reproduced. While Ap for the 
single domain lattice shows only a weak magnetic field 
dependence and stays always smaller than 5%, it shows 
peaks at the commensurability resonances and a broad 
shoulder extending from « 250 mT to « 500 mT, with a 
maximum up to 10 % in the GB lattice. Moreover, Ap 
has a minimum close to B — 200 mT, where the abso- 
lute values can even get negative. Fig. [IJb) shows that 
the temperature dependence of Ap between 2 K and 16 K 
is quite weak, suggesting a classical origin. Additional 



measurements at T = O.IK in a different cooldown (not 
shown) revealed no significant change of the effects as 
compared to 2 K. 

We proceed by interpreting our experimental results 
with the help of semiclassical simulations based on the 
Kubo formalismi^ The magnetoconductivity tensor was 
calculated in the lattice geometries as patterned in the 
heterostructure. The antidots were modeled as hard-wall 
circles with diameters similar to deff , namely 460 nm. 
For each value of the magnetic field, we computed the 
trajectories of 10^ electrons at the experimentally de- 
termined Fermi level of Ep = 10 meV in this potential 
landscape. The step width for the magnetic field was 
AB — 2mT. The electrons start at random directions 
of the Fermi velocity and at random positions within an 
area of 10 x 20 antidots. This large start area proved to 
be necessary to capture the transport properties of the 
whole GB within acceptable computation times. The 
electron trajectories are computed to a length of 75 pm 
from which we calculate the velocity correlation func- 
tion C^j{t,B) = {vi{t,B)vj{Q,B)), {i,j = x,y) for every 
value of Bi^ Here, the brackets denote averaging over 
all trajectories. The components aij{B) of the magne- 
toconductivity tensor for a degenerate two-dimensional 
electron gas follow from 
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{v,{t,B)v,{0,B))e-'/^"dt (2) 



Here, m* — 0.067me represents the effective electron 
mass in GaAs. In order to compare our experimental 
results with the simulations, the conductivity tensor is 
converted into a resistivity tensor. 

The simulations agree reasonably well with the main 
structure of the measured Ap{B), as can be seen in Fig. 
W[,c). The peak around 0.1 T, the minimum close to 0.2 T 
followed by a peak around 0.4 T are reproduced, with an 
amplitude close to the measured one. Since the simula- 
tions are carried out at zero temperature, this agreement 
supports our observation that thermal smearing is rela- 
tively unimportant. The simulated resistivity anisotropy 
for a single domain lattice, both without disorder as well 
as with isotropic disorder with a standard deviation of 
a = 0.1a (not shown), does not show any structure, in 
accordance with experimental results obtained by us as 
well as by others J^ However, the shift of the commen- 
surability resonance around one antidot towards larger 
magnetic fields as compared to the single domain lattice 
is not observed in the simulations. We speculate that 
in addition to the qualitative argument given above, a 
more realistic confinement potential is necessary to gen- 
erate this effect. 

We cannot resort to the conventional method of 
Poincare sectionsS, for further analysis, since our lattice 
containing the GBs is not periodic. We therefore follow 
a more qualitative approach and inspect the longitudi- 
nal velocity correlation functions of individual trajecto- 



ries, given by Cu{B) = / Wi(i, B)ui(0, B)e-*/^°di. For 
different magnetic fields, we identified trajectories tfiat 
contribute most to Ap. Sucli trajectories were found by 
inspection of the differences between Cyy{B) and Cxx{B) 
in the magnetic field intervals where the maximum of Ap 
occurs, namely around B = 400 mT. The result is exem- 
plified in Fig. [3l where we plot 300 trajectories starting 
from one unit cell. The GBs arc visualized by represent- 
ing the coordination numbers z of the antidots in gray 
scale. Coordination numbers of 5 and 7 are present in 
the GBs. The diffusion cloud has an anisotropic shape, 
in contrast to the corresponding cloud for a single domain 
lattice (not shown). Rather, the GBs form weak diffu- 
sion barriers, and characteristic branches of trajectories 
are visible that extend along the GBs. 




FIG. 3: Main figure: the diffusion cloud formed by 300 
electron trajectories that started from one unit cell (the 
white area at the center of the array) at a magnetic field 
of B = 400 mT. For better visibility, only every SO*** trajec- 
tory point, corresponding to a trajectory length of 250 nm, is 
shown. Inset: zoom of one trajectory bound to a GB. Note its 
localization at antidots with a coordination number different 
from 6 (the gray tone of an antidot represents its coordination 
number. Dark gray: z = 7, light gray: z = 6, and medium 
gray: 2 = 5). The dashed circle indicates the edge of an 
isotropic diffusion cloud. 

It is found that most of the trajectories in the branches 
are localized at a GB for a long time. The inset in Fig. 
13] shows one example of such a trajectory. It can be seen 
that the electrons show a tendency to circle around anti- 
dots with coordination numbers z of 5 or 7 for some time, 
before they move on to a neighboring antidot, which of- 
ten has a z 7^ 6 as well. For the simulations with the 
parameters of our samples, this behaviour is typical for 
magnetic fields above 0.3 T. At smaller magnetic fields, 
Ap is significantly smaller and even can become negative 
around B « 0.2 T, which indicates that the electrons dif- 
fuse better across the GB than along it. However, we 



could not relate this behavior unambiguously to a par- 
ticular type of trajectory in the corresponding diffusion 
clouds, probably because the effect is too weak. 



IV. SUMMARY AND CONCLUSION 

Longitudinal resistances parallel and perpendicular to 
artificially generated grain boundaries in a hexagonal an- 
tidot lattice imposed on a 2DEG in a semiconductor het- 
erostructurc have been measured. We observe that in 
comparison to a single domain lattice, the commensu- 
rability resonances broaden and develop a pronounced 
anisotropy. This suggests that the disorder across the 
grain boundaries is stronger than along them. The ef- 
fect of the grain boundaries on the electron trajectories 
has been visualized by numerical studies and consists of a 
slight localization of the electrons along the grain bound- 
aries for most, but not all, magnetic fields. In a small 
magnetic field interval in between the resonances around 
one antidot and seven antidots, both the experiments and 
the simulations suggest that an enhanced diffusion across 
the grain boudaries is possible as well. 

We simulated various GB types and could not iden- 
tify characteristic microconfigurations of antidots at the 
boundary. Therefore, we believe that the effects re- 
ported are qualitatively independent of the GB type. Fu- 
ture work could address several issues in further detail. 
Our choice of the ratio antidot diameter/lattice constant 
was based on a heuristic argument in combination with 
our experimental constraints. It is unclear how the ef- 
fects change with this parameter. Also, the influence of 
the modulation strength is an open question. In par- 
ticular, it will be worth looking at weakly modulated 
systems, which are more likely to be obtained by self- 
organizational schemes J^ The exact shape of the antidot 
potential, however, seems of minor relevance, since our 
samples have large lateral depletion lengths and with it a 
soft potential, and nevertheless the measured resistance 
anisotropics agree well with the simulations carried out 
on hard- wall arrays. 

In superlattice structures prepared by self-organization 
where GBs might matter, their type, density and orien- 
tation is per se at best marginally under control, and in 
samples containing many grains, the observed broaden- 
ing of the resonances will be isotropic. However, the dis- 
order generated by the GBs is non-Gaussian and locally 
anisotropic and thus different from random disorder. Our 
results indicate that even a high density of grain bound- 
aries does not destroy the specific transport properties of 
the superlattice, at least in the classical regime. 
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